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In this paper, we analyze the physical properties of dark mottles in the chromospheric network using two dimensional 
spectroscopic observations in Ha obtained with the Gottingen Fabry-Perot Spectrometer in the Vacuum Tower Telescope 
at the Observatory del Teide, Tenerife. Cloud modeling was applied to measure the mottles' optical thickness, source 
function, Doppler width, and line of sight velocity. Using these measurements, the number density of hydrogen atoms 
in levels 1 and 2, total particle density, electron density, temperature, gas pressure, and mass density parameters were 
determined with the method of Tsiropoula & Schmieder (1997). We also analyzed the temporal behaviour of a mottle 
using cloud parameters. Our result shows that it is dominated by 3 minute signals in source function, and 5 minutes or 
more in velocity. 
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1 Introduction 

The main fine structures of the quiet chromosphere emanat- 
ing from the network are called mottles when seen on the 
disk and spicules when seen at the limb. Mottles appear as 
dark or bright structures, which together form two kinds of 
groups, namely rosette and chain. The mottles in a rosette 
spread radially outwards from a bright center, while in a 
chain they all point in the same direction. The dark mottles 
are thin and elongated structures, whereas the bright mot- 
tles are small and roundish, lying at a lower height than 
the dark ones. These structures presumably extend along the 
magnetic fields and cover between 7" and 10" with widths 
smaller than 1" and have lifetimes of the order of 10 min or 
more. The spatial relationship between the bright and dark 
mottles is given in detail by Zachariadis et al. (2001) (see 
also references therein). 

There are a number of studies on the physical properties 
of chromospheric fine structures using spectrometric obser- 
vations as well as various theoretical models and simula- 
tions (Sterling 2000; Tziotziou et al. 2003; Al et al. 2004; 
De Pontieu et al. 2004; Rouppe van der Voort et al. 2007). 
Assuming an optically thin absorbing cloud in front of a ra- 
diation source, the standard cloud model (Beckers 1964) has 
been extensively used to extract intrinsic line formation pa- 
rameters such as source function, optical thickness, Doppler 
width and line of sight velocity (see review by Tziotziou 
2007). 

In this paper, we analyzed chromospheric features ob- 
served in Ha line with very good spectral profile sampling 
to derive the features' physical properties in terms of stan- 
dard cloud modeling. Using the method of Tsiropoula & 
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Schmieder (1997), additional physical parameters were also 
derived from the inferred cloud model parameters. We also 
addressed the dynamics of a dark mottle by the analysis of 
cloud velocity and source function variations. Our goal is 
to define the physical conditions inside dark mottles and 
to investigate the oscillatory behaviour by taking advantage 
of high resolution imaging spectroscopic capabilities of the 
Gottingen Fabry-Perot spectrometer and comparing our re- 
sults with previous studies. 

2 Observations and Data Analysis 

The data were obtained with the Vacuum Tower Telescope 
(VTT)'s Gottingen Fabry-Perot spectrometer, which is based| 
on two Fabry-Perot Interferometers (Koschinsky et al. 2001). | 
In 2002, a short time series of 60 wavelength scans of the 
Ha line were taken from a network region near the disk cen- 
ter of the sun. A spacing of 125 mA between adjacent wave- 
length positions was chosen for the spectral scans. With the 
narrow -band channel, 8 images at each of 18 wavelength 
positions were obtained, which provides a wavelength cov- 
erage of approximately 2. 125 A. The exposure time was 30 ms| 
and the time interval between the start of two subsequent 
spectral scans was 49 s. The observed field of view of the 
raw data was 38". 4 x 28". 6, with a spatial scale of 0".l per 
pixel. During the observations obtained under good seeing 
conditions, broad-band images were taken simultaneously 
with narrow -band images. Our best scan reached a Fried pa- 
rameter of ro=21 cm. Dark scans, flat fields and scans with 
the continuum source were also taken for the data reduction. 

Broad-band images were restored by using the spectral 
ratio (von der Liihe 1984) and the speckle masking tech- 
nique (Weigelt 1977). Narrow-band images were reconstructed! 
by using Keller and von der Liihe's (1992) method. 
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Fig. 1 A sample of speckle reconstructed narrow-band scan and the broad-band image (lower right) belonging to the same 
field of view. The size of each image is 32" x 23". 6. The wavelength positions decrease in each image by 125 mA starting 
from 6563.8 A at the upper-left corner. The last image in the second row is the closest to the center of the line profile. 



As a result of this method, we obtained one reconstructed! 
narrow-band image for each wavelength position of the Ha 
line. Figure [T] displays the reconstructed narrow-band im- 
ages through one scan. The observed region contains almost 
parallel elongated dark mottles, which form the so-called 
'chains of mottles'. 

For each pixel in the field of view, Ha line profiles were 
constructed from intensity values of narrow band images 
obtained at 18 wavelength positions. 

3 The Cloud Model 

The shapes and amplitudes of stellar spectral line profiles 
reflect the physical properties of the regions they are formed 
in. Therefore, physical parameters such as chemical abun- 
dance, density, temperature, velocity, magnetic field, micro- 
turbulence, etc. of these regions can be inferred through the 
analysis of observed spectral line profiles. 

The standard cloud model (Beckers, 1964) has been used( 
extensively to determine the physical properties of solar chro-| 
mospheric fine structures. The model considers a chromo- 
spheric feature like a cloud above a uniform atmosphere 



described by a reference background profile. The approach 
works well for optically thin structures and gives estimates 
of source function S, optical thickness at the line center tq, 
Doppler width AAd, and line-of-sight velocity «los f° r the 
observed cloud. These parameters are assumed to be con- 
stant within the cloud along the line of sight. The observed 
contrast profiles are matched with the theoretical contrast 
profiles using the following formula 
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where Irj(A) is the reference profile emitted by the back- 
ground and r(A) is the optical thickness. The wavelength 
dependence of optical thickness is given by 



t(A) = T exp 



A - A e (l - vlos/c) 

AAn 



(2) 



where A c is the line center wavelength and c is the speed of 
light. The standard cloud model considers chromospheric 
fine structures to be completely isolated from the surround- 
ing atmosphere and illuminated from the underlying atmo- 
sphere. 
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Fig. 2 Distribution of cloud parameters. Negative values 
indicate a downward motion in the line-of-sight velocity. 

The result of the cloud modeling is based on fitting the 
observed contrast profile with the model from an iterative 
least square procedure for non-linear functions. The line 
core intensity and the line core position of the observed pro- 
file were taken respectively as initial values of S and A c for 
the iteration procedure. For tq and AAd the values 1 and 
0.3 A were used. 

In order to determine the physical parameters of dark 
mottles using the cloud model, we selected dark features 
in the field of view by masking pixels greater than 0.9 of 
mean value of the line center intensity image for each spec- 
tral scan. Then we constructed line profiles for each selected 
pixel and calibrated them relative to the quiet sun contin- 
uum. Finally, we calculated a contrast profile for each se- 
lected pixel using the mean of all profiles in the field of view 
as the background profile and applied the cloud model. 

4 Results and Discussion 
4.1 Results of Cloud Modelling 

Figure [2] displays histograms of cloud parameters inferred 
for all the time series for dark mottles. The source function 
distribution shows a peak near 0.11 (in units of the contin- 
uum intensity I c ), while the Doppler width's peak is close 
to 0.46 A. The peak of the optical thickness distribution is 
around 0.8, which indicates that dark mottles are mostly op- 
tically thin structures. The distribution of line-of-sight ve- 
locities varies between —30 and 30kms _1 , and is almost 
symmetric, which implies the presence of both downward 
and upward motions. It peaks around — 1.25kms _1 (down- 
ward). 

Figure[3]shows maps of the cloud model parameters ob- 
tained for one scan. The map of source function has mini- 
mum values at the center regions of individual dark mottles. 
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Fig. 3 Maps of the source function, Doppler width, opti- 
cal thickness at the line center and line of sight velocity. In 
the medium grey areas either the cloud model fit was not ap- 
plied or did not give meaningful results at all. In the Doppler 
shift map, the bright and dark patches indicate upflows and 
downflows, respectively. 



Through the edge of the structures, the source function val- 
ues are close to the line center intensity of the background 
profile. We observe that regions of high Doppler width val- 
ues have low optical thickness and vice versa. This anticor- 
relation was also noticed by Alissandrakis et al. (1990) and 
Tsiropoula et al. (1993). On the map of the Doppler shift, 
downflows and upflows are always present along mottles 
while material is mostly descending near their footpoints. 

Table Q] shows values of the cloud parameters reported 
by various authors. It can be seen from the table that our re- 
sults, inferred for all the time series are in good agreement 
with previous measurements. Several intrinsic parameters 
can be determined by using the cloud model parameters. 
Assuming 15 km s -1 for microturbulent velocity, £ t , we can 
calculate temperature from the deduced Doppler width val- 
ues; 



AAn = — 



An 2kT 



m H 



(3) 



The relation between the optical thickness at the line 
center, the Doppler width and the number density in the sec- 
ond hydrogen level may be given: 



Ti'-e 2 /A 2 N 2 L , 
to = 1 — — a 



m e c c 



(4) 



inserting the constants in the equation above, we obtain: 



t = 1.38 x 10 2 



Nod 



AA 



D 



(5) 



The thickness of mottles estimated at the images was 
assumed to «350 km, considering a cylindrical structure, 
and the inclination was taken as «30° given for spicules 
by Heristchi and Mouradian (1992). Therefore the geomet- 
rical width, d of the structures under investigation, is equal 
to «700 km. Then using the same method as Tsiropoula 
& Schmieder (1997), we also obtained the electron density 
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N c , the total particle density Nn, the gas pressure p, the to- 
tal column mass M, the mass density p, and the degree of 
hydrogen ionization xh- The relations can be expressed as: 

N c = 3.2 lO 8 ^/^ 
7V H = 5 10 8 lO ' 5 losN2 
p = k(N c + 1.0851 N H ) T 

M = (N H m H + 0.0851 N H x 3.97m H ) d 

M 
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We assumed histograms of the parameters to be nor- 
mally distributed, so we fit them with a Gaussian curve sep- 
arately. We took peak values of Gaussian distributions as a 
mean value and their la widths as a standard error. In Ta- 
ble|2] we present the results obtained from all the time series 
and compare these with the values found by Tsiropoula & 
Schmieder (1997) and Tsiropoula & Tziotziou (2004) for 
dark mottles. The results are generally in agreement with 
previous computations, although some variations do exist. 
These differences may arise from various factors. First, the 
atmospheric seeing when the observation was performed, 
has a non-negligible effect on the determination of cloud 
parameters, especially on the Dopplerwidth (Tziotziou et 
al., 2007). Therefore, it causes an uncertainty in the calcu- 
lation of the pressure and temperature, which depends on 
the Dopplerwidth and the microturbulent velocity, which is 
assumed. Second, any uncertainties raised at the determina- 
tion of the geometrical thickness of the structures, d, will 
be propagated to the uncertainties in the calculations of N2, 
N e and M. Finally, differences in the inferred values of To 
can cause variations in these parameters. 

We calculated the pressure inside the dark mottles as 
0.32dyncm~ 2 . This result is in agreement with Heinzel & 
Schimieder (1994), who concluded that classical cloud model] 
can be applied to low pressure structures (< 0.5 dyncm -2 ), 
assuming a constant source function and a rather low opac- 
ity. 

Previous studies of spicules show that they have tem- 
peratures of about 5 000- 15 000 K, and densities of about 
3 x 10 -13 g cm' 3 . The physical properties we inferred from 
cloud modeling here for the mottles are in good agreement 
with the values found earlier for the spicules showing fur- 
ther evidence that these two features have significant simi- 
larities (Table|2]l. 

4.2 Temporal Analysis of a Mottle 

In order to study the oscillatory behaviour of chromospheric 
features, we performed a Fourier analysis to the fluctua- 
tions of cloud model's source function and Dopplershift for 
a long horizontal mottle's region. This region which is cen- 
tered at (8", 8") is shown with a white rectangle in Figured 
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Fig. 4 The mean image of Ha line center intensity series 
(top). White rectangle inside the image marks the mottle se- 
lected for the Fourier analysis. The distribution of optical 
thickness for the mottle subjected to timing analysis, (bot- 
tom). 



We selected this region because it could be well isolated 
from other structures and the cloud model fit it well for all 
of the time series. 

Before computing the power spectra, we averaged the 
signal in the selected mottle's region and then applied a 
polynomial trend removal of the 4th order to the time se- 
ries of source function and velocity. Consequently, we cut 
down the power at frequencies lower than 0.7 mHz, which 
corresponds to half of the total duration of our observation. 
This procedure also affected the amplitude of the power at 
the frequencies lower than 1.3 mHz, but the power at high 
frequencies remained unaffected. 

In Figure|5] we show the power spectra of Cloud Doppler-| 
shift and source function variations in the considered mot- 
tle's region. Power distribution of source function peaks at 
high frequencies between 5 and 8 mHz with a maximum 
around 5.7 mHz (i.e periods of 3 minutes) which corresponds! 
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Table 1 Cloud model parameters for dark mottles given by different authors. 



Authors 


S 


v (km s i ) 




AX D (A) 


Beckers 1968 


123 


- 


1.4 


0.5 


Bray 1973 


130 to 160 


-9 to 7 


1.0 


0.5 


Grossmann-Doerth & 


130 


-8 to 8 


1.1 


0.45 


von Uexkull 1977 










Tsiropoula et al. 1993 


163.3±14.3 


-0.26±6.6 


1.8±1.1 


0.37±0.1 


Lee et al. 2000 


0.16 I c 


-2.8 


2.2±0.5 


0.55 


Tziotziou et al. 2003 


0.15 I c 


-0.1 


0.9 


0.35 


Al et al. 2004 


0.14 I c 


0.18 


1.58 


0.44 


This work 


0.11 I c 


-1.25 


0.8 


0.46 



Table 2 Pyhsical parameters for dark mottles given by different authors. 



Parameters 


Tsiropoula & Schmieder 


Tsiropoula & Tziotziou 


This work 




(1997) 


(2004) 




iVi(cm- 3 ) 


(1.6±0.8)10 iu 




(3.0±0.7)10 1U 


N 2 {cm-' i ) 


(1.4 ± 1.1)10 4 


(4.2 ± 2.0)10 4 


(4.2 ± 1.9)10 4 


N e (cm- :i ) 


(3.4 ± 1.5)10 10 


(6.4 ± 1.6)10 10 


(7.0±1.6)10 10 


N H {cm~ 3 ) 


(5.1 ± 2.1)10 10 


(9.9±2.5)10 10 


(10.8±2.5)10 10 


T{K) 


(1.0±0.8)10 4 


(1.0±0.3)10 4 


(1.2±0.9)10 4 


P(dyncm 2 ) 


0.15 ±0.1 


0.24 ±0.1 


0.32 ±0.2 


M (gr cm~ 2 ) 


(2.2±0.4)10~ 5 


(2.2±0.6)10 -5 


(1.6±0.4)10" 5 


p (gr cm -3 ) 


(1.1 ±0.5)10" 13 


(2.2±0.6)10~ 13 


(2.4±0.6)10~ 13 


XH 


0.65±0.1 


0.65 


0.64 ±0.1 



to the acoustic cut-off frequency. The second significant powerj 
with a lower amplitude is around 7.2 mHz. De Pontieu et 
al. (2007) also reported a significant 3 minute periodicity 
in the intensity variations of the overlying loops emanat- 
ing from a network. They considered that the lower opacity 
of these structures allows glimpses of internetwork dynam- 
ics and oscillations underneath. In order to check this sce- 
nario, we looked into the distribution of optical thickness 
for the considered mottle (Figure |4j. The peak is around 
0.6, which is smaller than the peak of general distribution 
of optical thickness presented in Figure [2] This shows that 
the long mottle under consideration has a lower opacity and 
also confirms the suggestion of De Pontieu et al. (2007). In 
the power spectrum of source function, significantly lower 
powers are seen at low frequencies around 2.2 and 3.5 mHz. 

The cloud Dopplershift power spectrum shows three well| 
separated peaks at lower frequencies. The most significant 
peak lies around 2.2 mHz, i.e. at periods of 7 minutes. The 
other two are around 1.3 mHz and 3.5 mHz (i.e periods of 
5 minutes). We should note that magnitude of the power 
around 1.3 mHz is strongly affected by the procedure of 
trend removal, so we can not form a conclusion about its 
relevance. Also, peaks in the 3 minute range are visible but 
they have remarkably less significance. 

Similar to our results for the periodicities in velocity 
variations, Tziotziou et al. (2004) found a dominant period 
in the 5 minutes (3.3 mHz) for intensity and velocity varia- 
tion of dark mottles. De Pontieu et al. (2007) reported that 
the mottles are dominated by oscillatory behavior with peri- 
ods around 5 minutes and longer. More recently, Tsiropoula 



ft al. (2009) found a significant peak around the 5 minute 
range for the mottles' region. 

This strong peak around 5 minutes which was observed 
in the Doppler velocity power spectrum has long been at- 
tributed to p-mode oscillations generated at the photosphere 
propagate upward in and around magnetic flux concentra- 
tions of a network (Giovanelli et al. 1978). Several stud- 
ies based on simulations and observations have suggested 
that low frequency (<5 mHz) magnetoacoustic waves might 
propagate into the chromosphere along magnetic field lines 
which are significantly inclined with respect to the surface 
gravity (Bel & Leroy 1977; Suematsu 1990; De Pontieu et 
al. 2004; Jefferies et al. 2006). In order to study the propa- 
gation characteristics of waves at different heights in atmo- 
sphere, it is essential to compute phase difference and coher- 
ence spectra. In a future paper, we will concentrate on this 
issue in detail using the lambdameter method (Tsiropoula et 
al., 1993; Al et al., 2004) to compute intensity and Doppler- 
shift image series for different widths through the Ha pro- 
file. 

In this work, using high spectral and spatial resolution 
spectro-imaging observations, we constructed Ha line pro- 
files for every pixel of mottles in the field of view. We fit the 
resulting profiles with the cloud model and infered physi- 
cal parameters of the dark mottles, which agree with previ- 
ous studies. We also performed a Fourier analysis on fluc- 
tuations of the cloud velocity and source function to detect 
the oscillatory behaviour of the mottles. However, we could 
only concentrate on one mottle in the field of view since the 
cloud model can not fit the all of the profiles. Our analysis 
shows that the mottle is dominated by 3 minute periods (the 
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Fig. 5 Power spectra for the cloud source function and 
Dopplershift fluctuations of the mottle {top and bottom, re- 
spectively). 
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